NIT-221 
NT014 3US 



Title of the Invention 

A METHOD FOR MEASURING DISTANCE AND POSITION USING SPREAD 
SPECTRUM SIGNAL, AND AN EQUIPMENT USING THE METHOD 



Inventors 



Katsuhiko TSUNEHARA 
Nobukazu DOI 
Mikio KUWAHARA 
Tomoaki I SHI FUJI 



- 1 - 



TITLE OF THE INVENTION 

A METHOD FOR MEASURING DISTANCE AND POSITION USING 
SPREAD SPECTRUM SIGNAL, AND AN EQUIPMENT USING THE METHOD 

5 BACKGROUND OF THE INVENTION 

The present invention relates to terminal equipment 
for measuring its own position, particularly to equipment 
for measuring distances and positions using the radio waves 
emitted from base stations fixed on the ground, including 

10 CDMA base stations. 

The principles of distance measurement using a spread 
spectrum signal are described using Fig. 9. The station 
for transmitting the spread spectrum signal transmits this 
signal in send timing 400. The aforementioned receiving 

15 station receives the spread spectrum signal and obtains 
receive timing 401 . Differential time 402 between receive 
timing 401 and send timing 400 is detected as the 
propagation time of the spread spectrum signal. The 
distance between the transmitting station and the 

20 receiving station can be calculated by multiplying 

differential time 402 by the velocity of light. Because 
of the principles described above, distance measurement 
using a spread spectrum signal requires the measurement 
of receive timing 401 at the receiving station. 

25 Next, the principles of position measurement using a 
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spread spectrum signal are described. The distances to 
individual transmitting stations are measured by the 
receiving station, subject to the principles described 
above. The use of the thus -obtained distances between the 
5 receiving station and each base station and of the 

positions of the base stations enables the position of the 
receiving station to be detected by solving the equation 
where the position thereof is taken as an unknown quantity. 
Details of one such detection method are disclosed in, for 

10 example, Japanese Laid-Open Patent Publication No. Hei 
7-181242 (1995). 

To use spread spectrum signals for conducting distance 
or position measurements in this way, it is necessary to 
measure the receive timing of the aforementioned spread 

15 spectrum signal at the terminal equipment. In Japanese 
Laid-Open Patent Publication No . Hei 7-181242 (1995), the 
following method for measuring such receive timing is 
disclosed: the correlation values between the received 
signal and the predetermined code series for creating 

20 spread spectrum signals (hereinafter, collectively called 
the PN code) are calculated for each receiving event, and 
a profile is created that shows the values corresponding 
to the correlation values in each receiving event 
(hereinafter, this profile is called the delay profile); 

25 wherein an epitomized diagram of the delay profile is shown 
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as 1 in Fig. 10, and the timing where the correlation value 
becomes a maximum in the delay profile is searched for and 
the corresponding timing is detected as the timing in which 
the spread spectrum signal is received. In the example of 
5 Fig. 10, "t prev " is the receive timing. 

Summary of the invention 

During distance measurement and position measurement, 
it is important to measure the receive timing of the signal 

10 wave that has first arrived at the terminal equipment, 
namely, the incoming wave of the minimum propagation delay 
time. Consider the case that as shown in Fig. 11, a 
plurality of spread spectrum signals from a single spread 
spectrum signal transmitting station are passed along 

15 different propagation routes and received at terminal 
equipment as incoming waves 1 and 2 different in both 
propagation delay time and signal intensity. In this case, 
the delay profile received takes the shape of delay profile 
12, a combination of delay profiles 10 and 11 corresponding 

20 to incoming waves 1 and 2, respectively. In this case, only 
receive timing 22 of incoming wave 2 can be detected with 
the prior art. In the example of Fig. 11, since incoming 
wave 1 has the minimum propagation delay time and is 
received in timing 21, receive timing for the incoming wave 

25 of the minimum propagation delay time cannot be measured 
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using the prior art. As a result, receive timing 
measurement errors occur and this makes accurate distance 
or position measurement impossible. 

For these reasons, the use of the present invention 
5 enables the distance between a signal transmitting station 
and a signal receiving station to be measured by creating 
a delay profile from the signal wave received from the 
signal transmitting station, then taking the startup 
timing of the delay profile as reference timing, and 
10 detecting the timing delayed by a predetermined value 
behind the reference timing. 

To measure position, it is necessary to calculate the 
foregoing reference timing for at least three signal 
transmitting stations, then calculate the differences in 
15 send timing between the corresponding signal transmitting 
stations, and detect the position of the signal receiving 
station from the respective relative time differences. 

BRIEF DESCRIPTION OF DRAWINGS 

20 Figure 1 is a structural diagram of terminal equipment, 

the first embodiment of the present invention; Figure 2 
is a flowchart of the receive timing measurement algorithm 
used in the present invention; Figure 3 is a structural 
diagram of the delay profile creating section; Figure 4 

25 shows an example of a delay profile; Figure 5 shows the 
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first structural example of the first threshold value 
calculation section; Figure 6 shows the second structural 
example of the first threshold value calculation section ; 
Figure 7 shows the first structural example of the second 
5 threshold value calculation section; Figure 8 shows the 
second structural example of the second threshold value 
calculation section; Figure 9 is a diagram explaining the 
principles of distance measurement; Figure 10 is an 
epitomized diagram of a delay profile; Figure 11 is an 
10 epitomized diagram of the delay profiles created when two 
incoming waves are present. 

DETAILED DESCRIPTION OF THE INVENTION 

The receive timing measurement algorithm used in the 
15 present invention is described using the flowchart shown 
in Fig. 2, and an example of the delay profile shown in 
Fig. 4. 

In first step 500, the correlation value between the 
received wave and the PN code is calculated and delay 
20 profile 202 is created. 

In step 501, threshold value 206 required for making 
a distinction between incoming waves and noise 
(hereinafter, this threshold value is called the first 
threshold value) is calculated in delay profile 202. At 
25 this time, if in delay profile 202, the correlation value 
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exceeds the first threshold value 206, this threshold value 
is used to judge that an incoming wave is present in the 
particular timing, and this threshold value is 
sufficiently greater than the noise level. 
5 In step 502, among all the timing that the correlation 

value becomes equal to the foregoing first threshold value 
206, only the earliest receive timing 205 is detected 
(hereinafter, the earliest receive timing is called the 
first threshold value timing). 

10 In step 503 , threshold value 207 required for detecting 

the timing in which the delay profile corresponding to the 
incoming wave is calculated (hereinafter, this threshold 
value is called the second threshold value) . At this time, 
the second threshold value 207 is used to detect the timing 

15 in which the delay profile is started up from the noise 
level, and this threshold value is practically equal to 
the noise level. 

In step 504, among all the timing that the correlation 
value becomes equal to the foregoing second threshold value 

20 207, only the receive timing 208 closest to and earlier 
than the first threshold value timing 205 is detected as 
reference timing. Reference timing 208, therefore, 
denotes the timing in which the delay profile corresponding 
to the incoming wave is started up from the noise level. 

25 In step 505, the timing 210 delayed by predetermined 
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value 209 behind the aforementioned reference timing 208 
is calculated as reference timing. This means that the 
incoming wave has arrived at the receiving station in 
receive timing 210. Theoretically, predetermined value 
5 209, under its noiseless state, has a tip value of 1.0. 
In actuality, however, since noise exists, an edge 
subsequent to the true leading edge is detected as rise 
timing. This timing difference should therefore be 
subtracted to obtain a value from about 0.7 to 1.0. 

10 During position measurement that uses spread spectrum 

signals, when this measuring method, as with one shown in 
Japanese Laid-Open Patent Publication No. Hei 7-181242 
(1995), is to be used to conduct measurements using the 
relative distance differences between each transmitting 

15 station and the receiving station, step 505 can be omitted 
and, instead, the reference timing 208 obtained in step 
504 can be taken as receive timing 210. 

The construction of the terminal equipment, one 
embodiment of the present invention, is shown in Fig. 1. 

20 The spread spectrum signal that has been received by 

antenna 100 is sent to signal receiving section 101, where 
the signal then undergoes high/medium- frequency receiving 
and baseband signal demodulation. The spread spectrum 
signal, after undergoing processing in signal receiving 

25 section 101, is further send to delay profile creating 
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section 102. The correlation value between the received 
spread spectrum signal and the PN code is calculated for 
each receiving event by delay profile creating section 102, 
which then creates a delay profile that shows the values 
5 corresponding to the correlation values in each receiving 
event . 

A structural example of delay profile creating section 
102 using a matched filter is shown in Fig. 3. In Fig. 3, 
matched filter 200 calculates the correlation value 

10 between the received spread spectrum signal and the PN code 
created by PN code generator 201, and sends to signal line 
110 the value corresponding to the correlation value. An 
example of a delay profile created by delay profile 
creating section 102 is shown as solid line 202 in Fig. 

15 4. In Fig. 4, horizontal axis 212 denotes receive timing 
and as the delay profile bring closer to the left of the 
horizontal axis, the receive timing becomes earlier, that 
is, the propagation delay time decreases. Vertical axis 
213 in Fig. 4 denotes correlation values. 

20 The delay profile that has been created by delay 

profile creating section 102 is then held in delay profile 
holding section 115 . Delay profile holding section 115 can 
be, for example, a memory. The delay profile, after being 
held in delay profile holding section 115, is sent to the 

25 first threshold value timing detection section 103, the 



first threshold value calculation section 105, reference 
timing calculation section 106, and the second threshold 
value calculation section 107. 

The first threshold value calculation section 105 
calculates the threshold value to be used for the first 
threshold value timing detection section 103. A 
structural example of the first threshold value 
calculation section 105 is shown in Fig. 5. In this figure, 
the maximum value searching section 300 sends the maximum 
correlation value (existing in receive timing 203) of the 
delay profile received via signal line 110 . Multiplier 320 
multiplies the maximum correlation value 310 and 
coefficient C 0 and sends the results to the first threshold 
value timing detection section 103 as the first threshold 
value 330 . Coefficient C 0 is set to about 0.1. This avoids 
the likely mis-recognition of a side lobe caused by the 
characteristics of the band limiting filter within signal 
receiving section 101 during the creation of a delay 
profile; the side lobe being equivalent to a maximum 
correlation value 310 of about 0 . 1 in terms of magnitude. 

Another structural example of the first threshold 
value calculation section 105 is shown in Fig. 6. In this 
figure, noise power estimating section 301 estimates noise 
power using the delay profile received via signal line 110, 
and generates an output of noise power 311. The following 



two methods are available to measure noise power: 

(1) Approximating all received signal power to noise 

power 

(2) Creating a profile repeatedly and calculating the 
5 dispersion in the peak correlation values of the profiles 

Method (2) above, although higher than method (2) in 
accuracy , requires a long measuring time . Method ( 1 ) above , 
therefore, is used in Fig* 6. 

Multiplier 320 multiplies the abovementioned noise 

10 power 311 and coefficient C x and sends the results to the 
first threshold value timing detection section 103 as the 
first threshold value 330. Coefficient C x is set to a value 
from about 10 to 100 for this reason: when the noise is 
considered to be white noise, momentary amplitude changes 

15 in accordance with the required distribution, and in this 
case, if the noise power is taken as the square of a, the 
probability where the momentary amplitude exceeds 3a is 
about 3/1000, which is sufficiently slow as the probability 
of an measuring error occurring, and thus since an 

20 amplitude of 3a is nine times the square of a in terms of 
power, C l needs only to be more than nine. 

In Fig. 6, output 116 of signal receiving section 101 
can likewise be used as the input of noise power estimating 
section 301. Also, the first threshold value calculation 

2 5 section 105 can have the structural components shown in 



both Figs. 5 and 6, and send the greatest of the threshold 
values calculated thereby, to the first threshold value 
timing detection section 103 as the first threshold value 
330. Or the first threshold value calculation section 105 
can have the structural components shown in both Figs. 5 
and 6, and send the smallest of the threshold values 
calculated thereby, to the first threshold value timing 
detection section 103 as the first threshold value 330. 

The first threshold value 3 30 received from the first 
threshold value calculation section 105 is used for the 
first threshold value timing detection section 103 to 
generate the earliest receive timing in which the 
correlation value becomes equal to the first threshold 
value 330. The operation of the first threshold value 
timing detection section 102 is described using Fig. 4. 
In Fig. 4, discontinuous line 206 represents the first 
threshold value 330 received from the first threshold value 
calculation section 105. The earliest receive timing 205 
that, in delay profile 202, the correlation value becomes 
equal to threshold value 206 is sent from the first 
threshold value timing detection section 103 to signal line 
111. 

The second threshold value calculation section 107 
calculates the threshold value to be used for reference 
timing calculation section 106. A structural example of 
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the second threshold value calculation section 107 is shown 
in Fig. 7. In this figure, the same components as those 
shown in Fig. 5 as the first structural example of the first 
threshold value calculation section 105, are each assigned 
5 the same number as that of each shown in Fig. 5 . Multiplier 
320 multiplies the maximum correlation value 310 sent from 
the maximum value searching section 300, and coefficient 
C 2 , and sends the results to reference timing calculation 
section 106 as the second threshold value 331 . Coefficient 

10 C 2 is set to about 0.1, which is based on data that was 
measured using an experimental machine. 

Another structural example of the second threshold 
value calculation section 107 is shown in Fig. 8. In this 
figure, the same components as those shown in Fig. 6 as 

15 the second structural example of the first threshold value 
calculation section 105, are each assigned the same number 
as that of each shown in Fig. 6. Multiplier 320 multiplies 
the noise power 311 sent from noise power estimating 
section 301, and coefficient C 3 , and sends the results to 

20 reference timing calculation section 106 as the second 
threshold value 331. Coefficient C 3 is set to about 7, 
which is based on data that was measured using an 
experimental machine. 

In Fig. 8 f output 116 of signal receiving section 101 

25 can likewise be used as the input of noise power estimating 
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section 301. Also, the second threshold value calculation 
section 107 can have the structural components shown in 
both Figs. 7 and 8, and send the greatest of the threshold 
values calculated thereby, to reference timing calculation 
5 section 106 as the second threshold value 331. Or the 
second threshold value calculation section 107 can have 
the structural components shown in both Figs. 7 and 8, and 
send the smallest of the threshold values calculated 
thereby, to reference timing calculation section 106 as 

10 the second threshold value 331. 

The second threshold value 331 received from the second 
threshold value calculation section 107, the receive 
timing detection results received from the first threshold 
value timing detection section 103, and the delay profile 

15 received from delay profile holding section 115 are used 
for reference timing calculation section 106 to calculate 
the reference timing for obtaining the receive timing of 
the incoming wave of the minimum propagation delay time. 
The operation of reference timing calculation section 106 

20 is described using Fig. 4. In Fig. 4, single-dot dashed 
line 207 represents the second threshold value 331 received 
from the second threshold value calculation sectipn 107. 
Reference timing calculation section 106 compares the 
correlation value and threshold value 207 in the receive 

25 timing 205 that has been received from the first threshold 
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value timing detection section 103. If both values 
mismatch, the receive timing is advanced and the 
correlation value and threshold value 207 in said receive 
timing are compared. This sequence is repeated until the 
5 correlation value and threshold value 207 have matched, 
and the corresponding receive timing is sent as an output. 
In the example of Fig. 4, receive timing 208 in which the 
correlation value and threshold value 207 match is sent 
as reference timing to signal line 112. 

10 The reference timing received from reference timing 

calculation section 106 via signal line 112 is used for 
receive timing calculation section 108 to calculate the 
receive timing for the signal wave that has first arrived 
at the terminal equipment, namely, the incoming wave of 

15 the minimum propagation delay time. The operation of 

receive timing calculation section 108 is described using 
Fig. 4. Timing 210 delayed by previously set timing 209 
behind the reference timing 208 that has been sent from 
reference timing calculation section 106 is detected as 

20 the receive timing for the wave of the minimum propagation 
delay time, and the detected receive timing is then sent 
to signal line 113. 

The above method when applied to delay profile 12 shown 
in Fig. 11 is described. The first threshold value timing 

25 detection section can send receive timing 24 by using the 
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appropriate first threshold value 330. Next, the 
reference timing calculation section can send receive 
timing 20 by using the appropriate second threshold value 
331, Furthermore, receive timing calculation section 108 
can detect receive timing 21 by first measuring beforehand, 
under an environment having only one incoming wave, timing 
difference 23 between all values from the startup timing 
of the delay profile to the maximum value thereof, and then 
using said timing difference 23 in receive timing 
calculation section 108. Receive timing 21 is the receive 
timing for incoming wave 1 , the signal wave that has first 
arrived. In other words, even if two incoming waves are 
received in overlapping form, it is possible to detect the 
receive timing for the signal wave that has first arrived. 

Based on the receive timing 113 sent from receive 
timing calculation section 108, calculations for distance 
measurement or position measurement are performed by 
distance/position measuring section 114. Distance/ 
position measuring section 114 can use, for example, the 
method disclosed in Japanese Laid-Open Patent Publication 
No. Hei 7-181242 (1995). 

During position measurement that uses spread spectrum 
signals, when this measuring method, as with one shown in 
Japanese Laid-Open Patent Publication No. Hei 7-181242 
( 1995), is to be used to conduct measurements using the 



relative distance differences between each transmitting 
station and the receiving station, processing by receive 
timing calculation section 108 can be omitted and, instead, 
output 112 of reference timing calculation section 106 can 
be connected to signal line 113 and the corresponding 
output value can be sent to distance/position measuring 
section 114. In this case, delay profiles are created 
using the signal waves received from at least three signal 
transmitting stations, and then the first and second 
threshold values are created for each such delay profile. 
Subsequently, the startup timing of each delay profile is 
detected and the differences in send timing between the 
corresponding signal transmitting stations are used for 
the receiving station to measure its position from the 
relative time differences between the signal transmitting 
stations . 

The present invention enables accurate detection of 
the receive timing for the first incoming wave arriving 
under the multi-path environment that a plurality of 
incoming waves are received in overlapping form. Thus, it 
is possible to minimize measurement errors at the terminal 
equipment that uses spread spectrum signals to conduct 
distance and position measurements. 



